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Abstract 

We investigate the effect of the current measurement of the neutral B s meson mass difference, 
AMb s , on SUGRA models which have non-zero values of the soft breaking terms {m\ L ##)23 and 
A23 at the GUT scale. We use non-zero values of these parameter to explain the B — ► Kir puzzle 
and find that even after satisfying the experimental result on AMg s and the branching ratio (BR) 
of b —> sj we still can explain the puzzle. Further we show that in this parameter space it is possible 
to accommodate the large BR of B — ► n'K and the current experimental data for CF asymmetries 
of B -> n'K and B -> <pK°. The predicted value of sin (2/3^)^0 is about 0.52 - 0.67. 
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Flavor changing b — > s transitions are particularly interesting for new physics (NP) 
searches using B meson decays. In the standard model (SM) these transitions can occur only 
at the loop-level so that they are particulary sensitive to NP effects. So far, a few possible 
indications to NP effects through b — > s transitions have been reported by experimental 
collaborations such as BaBar and Belle. Among them is the recent B — > Kit puzzle: i.e., 
discrepancies between the SM predictions and the experimental results for the direct and 
mixing- induced CP asymmetries and the branching ratios (BRs) in B — > Kir modes whose 
dominant quark level processes are b sqq (q = u,d) Q, 0, [J • The measurements of the CP 
asymmetries in Bd — > rj'K and B& — > (f)K modes as well as the rather large BR for B — > rj'K 
and B — > rjK also have drawn a lot of attention, due to their possible deviation from the 
SM predictions The (dominant) subprocess of these modes is the b — * sss transition. 

Recently, the CDF collaboration has reported a new result for another interesting ob- 
servable relevant to the b — > s transition: the mass difference between the neutral B s states 
that characterizes the B s — B s mixing phenomenon. The CDF result is \^ 

AM Ba = 17.33l£g(stat.) ± 0.07(syst.) ps" 1 . (1) 

The D0 collaboration has also recently provided a new result 

17 ps" 1 < AM Bb < 21 ps" 1 (90% C.L.) . (2) 



These experimental results are consistent with the SM estimation. Therefore, these new 
experimental results are expected to provide important constraints on NP beyond the SM 0. 
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dv these new results, some theoretical studies have been done to search for NP 
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In the SM, the mass difference in the B s system is given by 



am| m = 9M3l Mb ^ b b b jI \v tb v;f s ( Xt ) 



(3) 

2.463 Q- 



where the NLO short-distance QCD correction gives r) B = 0.552 and So(x t ) 
The non-perturbative quantities B Bs and f Bg are the bag parameter and the decay constant, 
respectively. The best fit for AM™ is given by HQ 



AMi M 



21.5 ± 2.6 ps" 1 [UTfit] , AMg 



SM 



21.71H 1- ; 



[CKMfitter] , (4) 



In a recent paper 



Q. 



this mass difference is found to be 23.4 ± 3.8 ps using HPQCD 



and JLQCD data for f Bs y B Ba . 



In this letter, we study the neutral B s meson mixing effect in supersymmetry (SUSY): 
specifically in the supergravity (SUGRA) model. Then, using the constraints obtained from 
AMb s , we focus on how to resolve all the possible current anomalies observed in hadronic 
B — > PP (P denotes a pseudoscalar meson) decays through the b — > s transitions, such as 
B — > Kir, B — > r/fT. The current experimental data is listed in Table 1. 

We consider the SUGRA model with the simplest possible non-universal soft terms which 
is the simplest extension of the minimal SUGRA (mSUGRA) model. In the SUGRA model, 
the superpotential and soft SUSY breaking terms at the grand unified theory (GUT) scale 



are given by 



W 

£soft 



Y U QH 2 U + Y D QH 1 D + Y L LH 1 E + /iH^, 
1 



J2 M a \ a K + B^iE x E 2 



+ (A U QH 2 U + A u QH l D + A LH 1 E) + H.c. 



(5) 



where E, U and D are respectively the lepton, up-quark and down-quark singlet superfields, 
L and Q are the SU(2)^ doublet lepton and quark superfields, and E\ 2 are the Higgs 
doublets, (pi and X a denote all the scalar fields and gaugino fields, respectively. 

The SUSY contributions appear at loop order. In our calculation, we do not use the 
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23j . We assume the 



m 



(QlL,U R r,D RR ) 



Sij + [A(Q LL ,U RR ,D RR ) 



mass insertion approximation, but rather do a complete calculation 
breakdown of the universality to accommodate the b — > s transition data. While we satisfy 
this data, we also have to be careful to satisfy other data, e.g., b — > sj. 
We use the following boundary conditions at the GUT scale: 

, tn , A^ = A (Y^ + AA^) , (6) 

where i,j = 1, 2, 3 are the generation indices. The SUSY parameters can have phases at the 
GUT scale: Mi = \M 1/2 \e idl (the gaugino masses for the £7(1), SU(2) and SU(3) groups, 
i = 1, 2, 3), A Q = \A \e ia A and /i = |/x|e i6,M . However, we can set one of the gaugino phases to 
zero and we choose 9 2 = 0. The electric dipole moments (EDMs) of the electron and neutron 
can now allow the existence of large phases in the theory 2J]. In our calculation, we use 0(1) 
phases but calculate the EDMs to make sure that current bounds (\d e \ < 1.2 x 10~ 27 ecm |25l| 
and \d n \ < 6.3 x 10~ 26 ecm j26j) are satisfied. 

We evaluate the squark masses and mixings at the weak scale by using the above boundary 
conditions at the GUT scale. The RGE evolution mixes the non-universality of type LR (A 



TABLE I: Experimental data on the CP-averaged branching ratios {B in units of 10 6 ), the direct 
CP asymmetries (Acp), and the effective sin(2/3) (/5 is the angle of the unitarity triangle) for 
B — > PP decays |l|. 



BR 


Average 


CP asymmetry 


Average 


B(B ± 




24.1 ± 1.3 


Acp{K u ?r 


—0.02 ± 0.04 


B(B ± 


i^vr ) 


12.1 ±0.8 




±0.04 ± 0.04 


B(B° 




18.9 ±0.7 


„4 C p(-^ T ) 


-0.115 ±0.018 


B(B° 


K vr ) 


11.5 ± 1.0 


A7P(-F^°) 


±0.001 ±0.155 








sin(2^ eff ) JCs7r o 


±0.34 ±0.29 


B(B ± 


-» r/if ± ) 


69.7±l;f 


sin(2/3 eff )^ K o 


±0.50 ±0.09 








sin(2/3 eff )^ 


±0.47 ±0.19 



terms) via drriQ 2 LL RR / dt oc A\^A u ^d) terms and creates new LL and RR contributions at 
the weak scale. We then evaluate the Wilson coefficients from all these new contributions. 
We have both chargino and gluino contributions arising due to the LL, LR, RR up type 
and down type squark mixing. These contributions affect the following Wilson coefficients 
(73— (79, (77 7 and C$ g . The chargino contributions affect mostly the electroweak penguins (C7 
and C9) and the dipole penguins, where as the gluino penguin has the largest contribution 
to the dipole penguins due to the presence of an enhancement factor nig /nib (The gluino 
contribution also affects the QCD penguins). We include all contributions in our calculation. 

For calculation of the relevant hadronic matrix elements, we adopt the QCD improved 
factorization. This approach allows us to include the possible non-factorizable contributions, 
such as vertex corrections, penguin corrections, hard spectator scattering contributions, and 
weak annihilation contributions 27 1. 

The neutral B meson mass difference involves gluino and chargino diagrams in SUSY j3| . 
In mSUGRA, with universal boundary condition, the chargino diagram has the dominant 
contribution. Once we introduce mixing in the (2,3)-sector of the m\ L RR or soft 
breaking terms, the mass difference gets enhanced and we get large contributions from the 
gluino diagrams. The B — > ttK puzzle can not be solved using just the mSUGRA boundary 
condition. In order to explain the B — > puzzle, we have noticed that the flavor violating 

n 

terms in the (2,3)-sector of the soft breaking masses are needed pj. 



4 



-0.095 



34 35 36 37 38 39 

AM B ,/AM Bj 



FIG. 1: AcpiK^it*) versus AMb s I AMs d in the SUGRA model. The parameters are described 
in the text. 

In order to investigate the effect of the neutral B s mixing on the h — > s transitions, we 
first try to fit the B — ► Kir data using A^f, {^ll rr)^ a ^ the GUT scale. The constraint 
from the BR of b — > 37 is also included. We vary m 1( / 2 in the range (350 — 500) GeV 
[corresponding to gluino mass of (1 — 1.5) TeV], A = —800 GeV, &-(q ll ,u rr ,d rr ) ~ — 0.3, 
AA^ 3 ' d) = - 0.3, m = 300 GeV and tan (3=40. The A's also have 0(1) phases. The 
magnitudes of A's get reduced at the weak scale compared to the GUT scale since the 
squark masses get a contribution from mi/2 in the RGEs. 

In Fig 1, we plot Acp{K ± tt t ) versus AM^/AM^, where AMg d is the mass difference 
between the neutral states. The experimental value for the AMg d is 0.507 ± 0.005 ps _1 . 
In the SM, 



AM S SM M Ba , V ts 



AMf M Bd 



V, 



td 
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(7) 



where £ = fesV ^ . In the plot, we used £ = 1.18 and the CKM phase 7 = 61.1° ± 4.5° 



We find that the 2a experimental range about the central value of the ratio AMb s / AMs d = 
34.66 rules out a lot of model points. In order to extract the valid points, we include the 
error of £ = 1.23 ± 0.06 Q (consistent with the value of £ = 1.21±g;8fl using the JLQCD 
and the HPQCD calculations in Ref. Q|), and calculate the BRs and the CP asymmetries 
of different B — > Kn modes. We also calculate the BR of B — > rj' K and sin(2/3 eff )^/^o as 
well as sin(2/3 eff )^o for the allowed model points. 

The recent experimental data for the CP-averaged BRs of B — > Kit may indicate a 
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FIG. 2: R r -R n versus R r in the SUGRA model. 
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FIG. 3: Acpi^n") versus Acp^tt^) in the SUGRA model 



possible deviation from the prediction of the SM: 



Rr 



2B{B d 



1.00 ±0.09 



Rn 



B(B° -> X ± 7T^) 

2£(5° -> iW) 



0.79 ±0.08 . 



(8) 



It has been claimed that within the SM, R c m 



Rr. 



29 



30j. But, the data show the pattern 



R c > Rn, which would indicate the enhancement of the electroweak (EW) penguin and/or 
the color-suppressed tree contributions In Fig. 2, we plot R c — R n versus R c and find 
that R c > R n can be satisfied. 

Also, in the conventional prediction of the SM, Acp{K ± h°) is expected to be almost the 
same as Acp(K n^). In particular, they would have the same sign. However, the current 
data show that A C p{K ± 'K ) differs by 3.5a from Acp{K ± tc*)- In Fig. 3, we plot Acpi^rr ) 
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FIG. 4: sin(2/3 eff ) Ks7r o versus Acp(K±if?) in the SUGRA model. 
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FIG. 5: sin(2/3 eff ) Vjft: o versus the BR of £± r/'K ± in the SUGRA model. 

versus Acp{K ± ti t ) and find that the signs can be different for the points allowed by the 
neutral B mixing data. 

The predicted sin(2/3 eff )# s7r o is shown in Fig. 4. We find that the minimum value is 0.7. 
The present experimental data still have large errors so that future results will confirm/rule 
out our model. 

The experimental BRs of B(B — > rj'K) are large compared to the conventional SM pre- 
dictions. In Fig. 5, we plot sm(2/3 eS ) v r K o versus the BR of B ± — > r]'K ± . These decay modes 
get SUSY contributions since we are using non-zero values of {jn 2 LLRB )23 and A 2 3 and the 
BR gets enhanced. The values of sm(2(3 eS ) rj r K o is allowed by the experimental value which 
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TABLE II: The CP-averaged branching ratios (B in units of 10 6 ), the direct CP asymmetries 
(Acp), and the effective sin(2/3) for m l / 2 = 450 GeV, mo = 300 GeV and tan/3 = 40, IA23L1J = 
0.48, \AA*\ = 0.1, \AAl 



23 



0.3. 



BR 


Average 


UP asymmetry 


Average 


B(B ± 




23.8 




-0.03 


B(B ± 


-> K^ ) 


11.1 




0.013 


B(B° 


-» ifV) 


19.6 




-0.10 


B(B° 




11.4 


Acp{K°ir ) 


-0.11 








sm(2p cS ) Ksn0 


+0.8 


B(B ± 




72 


sin(2/? efl %, K o 


+0.6 



has a smaller error than that of sin(2/3 efr )^o. The value of sin(2f3 eS ) ( f ) x is predicted to be 
around (0.55 - 0.70) and the BRs of B ± -> QK*- and 5° -> </>fsr° are around (7 - 9) x 10~ 6 
and (6.5 — 8.5) x 10 -6 respectively in our calculation and B(b — > 57) is (2 — 4.5) x 10~ 4 . We 
also find that B(B — ► r/X) is around 3 x 10 -6 . The CP asymmetries for B ± — > ^if^ and 
i? 1 * 1 — > r]^K ± are —0.1 to 0.1 and close to 0, respectively. The Arg[M(12)sJ is less than 5° 
for our model points. 

It is possible to obtain a fit for the experimental results even without using m\ L con- 
tribution at al. The nonzero values of parameters generate the dipole penguin and 
the (Z-mediated) electroweak penguin diagrams. As a representative example, we present 
the BRs and CP asymmetries for a specific model point in Table II to show that all these 
different experimental results can be explained by one model point using just A?fa . The pa- 
rameters of the model point are given at the GUT scale by mi/2 = 450 GeV, mo = 350 GeV, 



An 



-800 GeV, AA 



23 



0.1 e" 2 - 0i , AA% 3 



0.48 e , and we choose tan/3 = 40. We find 



that the BRs and CP asymmetries are all within one sigma of the experimental results except 
for sin(2/3 eff )x S 7r° which is about 1.6a away (this deviation is lowered when we include m 2 LL 
contribution). The QCD parameters for this fit are qa = 2 and <pA = 2.75, where pa and 4>a 



are defined by X A = Jq f = (l + p A e^ A ) ln^ |27j. The ratio AM B JAM Bd is 34.3 for this 
model point. The EDMs are following: \d e \ = 2.48 x 10~ 29 e cm and \d n \ = 8.6 x 10~ 28 e cm. 
The BR of b -> s 7 is 4.2 x 10~ 4 . 

The origin of the (m 2 LL #^2)23 terms are natural in the grand unifying models which explain 
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neutrino masses. For example, if right handed neutrinos exist, SU(5) might generate the 
term Y v 5N5h, where 5 has d\ and the lepton L doublet, N is the singlet right handed 
neutrinos and 5# contains the SM Higgs doublet (along with the colored Higgs fields. Y v 
has a flavor structure in order to explain the neutrino masses and bilarge-mixing angles. 
Now these couplings introduce flavor violation to the soft masses (d c and I) via the RGEs, 



(x m 2 Y u Yj. In this model the m?. g terms for i ^ j can be generated 



. These terms get 



introduced between the GUT scale and the string scale due to the RGEs. One expects these 
flavor violating terms also in the SO(10) type models |32]. The right handed neutrinos there 
belong to the fundamental 16 representation of SO(10) and produce these flavor violating 
terms in the soft masses. The flavor structures of the Dirac and Majorana coupling arise from 
the neutrino mixing matrix. The terms (for % ^ j) also get contributions from the flavor 
structure of Y v due to the quark-lepton unification. Similar flavor violating effects in the 
soft terms are also present in the Pati-Salam type models [33]. In this case, the quark-lepton 
unification can happen at the intermediate scale and the flavor violating Majorana coupling 
/'04,i,2'03,i,2^io,i,3 (^4,1,2 contains right handed neutrinos along with the right handed quarks 
and leptons, Aio,i,3 is the new Higgs field ) is responsible for right handed neutrino Majorana 
masses. Now the RGEs involving these couplings between the intermediate scale and the 
grand unifying scale can easily introduce flavor violating terms in the squark and the slepton 
masses. 

In conclusion, we find that the current experimental results on the neutral B s meson 
mass difference have introduced strict constraint on the the SUGRA parameter space for 
flavor mixing terms A23 and {m\ L ^2)23 in the soft SUSY breaking terms. These flavor 
violating soft breaking terms are natural in the grand unifying models. In order to explain 
the B — > Ktx puzzle, A23 and (m 2 LL RR )23 are needed. We show that it is still possible to 
explain the B — > Kir puzzle even after satisfying the new Tevatron result on AM# s . The 
model used here contains three complex nonuniversal soft breaking terms (A23,ll, AA23 ), 
though an acceptable fit can be obtained using just AA^. This allows us to calculate 
19 observables of the B system (9 observables in the B — > ttK modes, 4 observables in 
B — > (f)K modes, 5 observables in the B — > rj^K modes and B(b — ► sj)). The future results 
on sin(2/3 eff )^ s7r o and Acp{K ± n T ) are crucial to probe this model. Finally, we find that the 
large B(B — > rj'K) can be explained in this parameter space with sin(2/5 eff )^xo near the 
current experimental result which is 2a away from sin(2/3 efr ) j/^k ■ 
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